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Abstract 



We consider the production and decay of the \2 charmonium state in polar- 
ized and unpolarized pp collisions at RHIC in the framework of an effective 
theory based on the QCD multipole expansion. We find that the angular 
distribution in the decay of the produced charmonium, \2 — ► J/ip + 7, in 
the unpolarized case allows us to distinguish clearly between the color singlet 
and color octet production mechanisms. Once the production mechanism is 
known, the angular distribution in the polarized case can be used to measure 
the polarized gluon distribution in the proton , AG(x). 
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The production of heavy quarkonia recently has attracted much interest. Tevatron data 
IjH have clearly shown that the traditional color singlet model 0] fails to explain the mea- 
sured J / ij) and ip' production cross sections, and have stimulated the development of a more 
general approach based on the non-relativistic expansion in QCD ||. This approach sys- 
tematically classifies, in powers of heavy quark velocity, the contributions of intermediate 
color singlet and color octet quark-antiquark states to quarkonium production. The cor- 
responding long-distance color octet matrix elements unfortunately cannot be computed 
analytically at present, and have to be either extracted from the data or computed on the 
lattice. 

In this Letter, we use the ideas of QCD multipole expansion |4| to formulate a simple ef- 
fective approach to treat both color singlet and color octet intermediate states in quarkonium 
production. Specifically, we consider the production of \2 charmonium state in polarized 
pp collisions at small Pf. This process has been considered previously by several authors 
Here we show that the measurement of the angular distribution in the decay of the 
produced charmonium, \2 - ► J/ip + 7, can be used both to distinguish between color octet 
and singlet production mechanisms and to measure the polarized gluon distribution in the 
polarized nucleon. The knowledge of the production mechanism is a necessary pre-requisite 
for the extraction of AG(x, Q 2 ) from the spin asymmetries in charmonium production. It 
will also eliminate many uncertainties that presently complicate interpretations of the data 
on charmonium production in nuclear collisions f|. 

We begin with a summary of our results. The underlying mechanism at the parton level 
is gluon-gluon fusion. First we consider X2 production in unpolarized pp collisions. The \2 
may be produced directly - the color singlet mechanism - with amplitude A 1 . As pointed 
out in Ref. 0, direct production at low p± produces with helicity ±2 only. Alternatively, 
the X2 can be produced in a two step process involving a color octet intermediate state. First 
two gluons annihilate to form a cc pair in a color 8 state. As we argue below, the J PC = 1 
cc configuration dominates. This decays by a color El transition to the \2 with helicity ±1 
only. We parameterize the production and decay process by an amplitude A s . A 1 and A 8 
are independent of the kinematic variables x\, x 2 and 9. The angular distribution of the 
photon relative to the beam axis in the decay \2 ~^ J + 7 determines the ratio of -4 s to 
A\ 

^ cx g( Xl , Ml)g{x 2 , M 2 ) + * cos 2 e)A 1 +(^-\ cos 2 9) A 8 } (1) 

Once A s / A 1 has been measured in unpolarized pp collisions, the same process can be 
used to measure Ag(x, M 2 ) at a polarized pp collider like RHIC. The angular distribution 
of the spin asymmetry in \ 2 production is given by, 

da^ - dan Ag(x 1 ,M 2 )Ag(x 2 ,M 2 ) ^ \ + \ cos 2 9 - g{| - \ cos 2 9} 

do-V + do-U ~ g{x u MD g(x 2 ,M*) X \ + \ cos 2 9 + § {§ - \ cos 2 9} U 
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Equations (P and (0) are our basic results. After a brief discussion of the restrictions on 
their validity, we present a derivation. We assume that the \2 is produced by gluon-gluon 
fusion. Light quark fusion is known to become important only at xp > 0.4 [|10|| . We ignore 



the gluon transverse momentum. It is easy to show that intrinsic transverse momentum 
generates corrections of order k^/M^. The most important kx correction is the possibility 
to produce the \2 with A = 0, with amplitude ~ k T /M 2 . kx effects can be eliminated 
altogether experimentally at some cost in rate, by limiting the total transverse momentum 
of the photon and dilepton to a small value. 

We work to lowest non-trivial order in perturbative QCD. Hard QCD corrections to 
charmonium production by gluon-gluon fusion have been studied in Ref. [|ll|]. They can be 
ignored for X2S produced at low kr- Our results follow from three observations: 

1. Direct production of the \2 in the color singlet state occurs only in the helicity ±2 
states. 

2. Production of the \2 through an intermediate color octet state yields %2's predomi- 
nantly with helicity ±1. 

3. Each X2 helicity state decays into J/ip + 7 with a characteristic angular distribution. 

Although (2) is the only ingredient not either in the literature or obtainable by elementary 
means, for completeness we review all three. 

We describe the \2 state of charmonium by a composite tensor field t/" 7 , which is sym- 
metric, 7]^ = rf^ 1 , traceless, 77^ = 0, and transverse, = 0. It is convenient to define 
the X2 wave function H I1U (P,A) with fixed momentum P and helicity A: 

(P,A|^|0) = H» U (P,A) with H^ = H U>1 , H£ = 0, P^H^ = 0. (3) 

In the rest frame of the X2, the transversality condition implies that H 00 = H°i = iP° = 0. 
The wave function in the rest frame therefore reduces as expected to an irreducible 3-tensor 
of rank 2, with 5 independent components describing different spin projections. 

We must consider two couplings of the \2 state to gauge vector bosons. First, the 
electromagnetic decay coupling to the photon plus J/ip and second, the QCD production 
coupling to a [1 ] 8 (cc) state and a gluon. Both are electric dipole transitions with effective 
interaction Lagrangians, 

£ n = ^iiTr{*tD,}^ ~ Tr{*t Ai/ }^. (4) 

where ip^ (^/i) is an interpolating field for the J/ip ([1 ] 8 (cc)) and (A M ) is the photon 
(gluon) vector potential. A M and are matrices in the octet representation of SU(3) co i or , 
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and Tr denotes a trace over color indices. The effective Lagrangians eq. correspond to 
the lowest order in the multipole expansion and are valid when the momentum p of the 
gauge boson (photon or gluon) is small compared to the inverse size of charmonium i? -1 : 

p « Br 1 . 

From eq. (f|) it is straightforward to evaluate the angular distributions in the decay 
X2 —>■ J/ip + 7 when \2 is produced in different helicity states. It is convenient to perform 
the calculation in the Xi res t frame where the amplitude corresponding to eq. (f|) can be 
written as, 

Aoce k \X^(X,)H kl (A), (5) 
where e's are the polarization vectors of J/ip and 7. Using the completeness relations, 

X; e m (X i) )e k \\ i) ) = 6 km , ]T e"(A 7 )e't(A 7 ) = 5 nl - k n k l ; (6) 

for the polarization vectors, where k is the unit vector in the direction of the photon's 
momentum, the cross section can be written down as 

^oc J2 A^A(xH^ kl (A)(5 ln -k l k n )H kn (A). (7) 

The evaluation of eq. (0) simplifies if one parameterizes iiP in terms of the unit vectors 
v(/i) that describe the angular wavefunction of a state with J = 1 and Mj = fi, 

W k {K) = X(^V|2A)v^( M )v fc (/iO. (8) 

The vector v(/i), can be represented in a Cartesian coordinate system as 

v(±l) = T^(ei±te 2 ), v(0)=e 3 . (9) 

The phase convention in eq. (|9]) corresponds to the familiar Cartesian representation of vector 
spherical harmonics; the v(/z) are normalized to unity, |v(/i)| 2 = 1. The decomposition of a 
vector w with respect to this basis is 

w = cos 9 v(0) + -p sin 9 v(-l) - -= sin 9 v(l) (10) 
v 2 v 2 

It is convenient to choose 9 as the polar angle relative to the spin quantization axis (which 
we identify with the hadron beam axis); the polar angle <fi does not play any role in our case. 

The relations eq. (|6|fl|9|JT0D now make the evaluation of the angular dependence of the 
cross section (|7|) trivial. For different projections of the X2's spin, we find 
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O)oci + |co^, 

H^ ±1 (#)oc--icos 2 #, 
w 4 4 

W/ 2 '°(#)oc§-^cos 2 #. (11) 

U — 

These decay distributions are rather distinct, a fact which can be used for experimental 
determination of the helicity state of the \2- 

Let us now turn to the mechanism of \2 production. We are concerned with two com- 
peting mechanisms - "color singlet" and "color octet" production. In the first, two gluons 
annihilate directly to the \ 2 - a color 1 cc state with J PC = 2 ++ . In the second, two gluons 
annihilate first to a color 8 cc state which then decays to the \ 2 ■ 

The helicity dependence of the color singlet mechanism was studied in || , who concluded 
that the \2 is produced only in the helicity ±2 state. Helicity is allowed by conservation 
laws but does not occur. We agree with this result. Our calculation proceeds as follows. 
First we compute the amplitude, M., coupling a \2 at rest with helicity A to a non-relativistic 
quark antiquark pair parameterized by a relative momentum q and two component spinors, 
U and V. The momentum space wavefunction of the X2 , 0(9), need not be specified except 
that it falls rapidly with q = \q\, 

^(A,g,f/,y) = 0(g)^(l/il / il2A)(v(/i)-gl(v( / iO-FW) (12) 

mm' 



The a and q in eq. (|i~2"|) reflect the fact that the \2 state is a spin triplet p-wave excitation 
of the cc pair. 

The cc pair produced by gluon-gluon fusion have a distribution in spin and momentum 
determined by the elementary tree diagrams for gg — > cc. Because this process involves gluon 
momenta p of the order of the charm quark mass, p ~ M c , we are justified in describing 
the annihilation process to lowest non-trivial order in QCD. From the Feynman diagrams 
it is straightforward to compute the q and a dependence of the gg — > cc amplitude, Af, (to 
lowest non-trivial order in q) as a function of the gluon polarizations, 

J\f(q, el, e 2 , U, V) oc U ] e x ■ q e 2 • a + e 2 ■ q e% ■ a - e x ■ e 2 q ■ a + e x ■ e 2 q ■ h a ■ k x V, (13) 

where we have dropped terms zeroth order in q which vanish when averaged over the wave- 
function of eq. (0). e*i and e 2 are the polarization vectors for the incoming gluons with 
momenta k\ = —k 2 in the center of mass; k\ is a unit vector in the direction of k x . 

The helicity amplitudes for gg — > \2 are determined by combining eqs. ( p!2[) and (fL3D , 



summing over quark spins, and integrating over q. Since we are interested only in the relative 
amplitudes for different helicities we ignore the integral over the magnitude of q, 
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K{A, el, e 2 ) oc / ^ £ A<t(A, g, [/, \/)Af(g, el, e 2 , £7, 1/) (14) 

From this it is straightforward to see that 1Z vanishes when A = 0. Simple helicity con- 
servation and symmetry arguments then suffice to show that the amplitudes for gluons of 
opposite helicity to form a \2 with helicity ±2 are equal (up to an inconsequential phase). 
This completes our confirmation of point (1) in our list. 

The discussion of color octet production proves simpler than color singlet production, 
because the helicity amplitudes are all determined by symmetry considerations once the 
dynamical mechanism is clear. The first step is the annihilation of two gluons into a color 
8 cc-state. The (cc) 8 state can have a variety of J quantum numbers. It then decays to 
the \2 by gluon emission. If the relative momentum of the (cc) 8 pair q is large compared 
to the inverse radius of quarkonium R~~ l , the formation of a bound state is very unlikely. 
Formation of quarkonium therefore requires a small invariant mass of the pair, and in this 
case the momentum of the radiated gluon p is typically small compared to -R -1 , so the 
QCD multipole expansion in powers of pR should be meaningful and convergent. These 
considerations strongly favor electric dipole gluon emission for the decay (cc) 8 — * X2Q- We 
are forced to conclude that the parent (cc) 8 state has J PC = 1 . The only other state that 
could couple by El radiation would be 3 . The 3 is a <i-wave internal excitation of the 
cc pair as opposed to the 1 state, which is an s-wave. 

Having established that the parent color-8 state is a 1 , we can use C\ and C\\ of eq. (|j) 
to extract the relevant helicity amplitudes. The logic is simple: The [1 ] 8 has helicity states 
±1 and 0. Angular momentum conservation dictates that only the helicity-0 state of the 
[1 ] 8 can be produced at rest or low momentum in gg fusion. When the helicity-0 [1 ] 8 
state decays to the [2 ++ ] 1 X2, it can only yield the \2 in the helicity ±1 state. This completes 
our argument for point (2) in our list. 

Our principal results, eq. (JT|) and (§), follow directly from the helicity selection rules (1) 



and (2) and the angular distributions of eq. (11). The attractive feature of this analysis 
is that the reaction mechanism can be determined by the study of the \2 decay angular 
distribution in the unpolarized case (eq. (|l|)), after which the gluon distribution in the 
proton can be determined by analysis of the polarization asymmetry of the decay angular 
distribution (eq. (fj)). 
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